function [P]=plate_pair(region,ws); 
%Identification of the MORVEL plate pairings for mid-ocean-ridge segments
%analyzed in this study. 

R(1).name='Pacific'; 
R(2).name='Cocos';
R(3).name='Nazca';
R(4).name='Antarctic';
R(5).name='Australian';
R(6).name='Juan de Fuca';

%late index/number 
R(1).index=0;
R(2).index=6;
R(3).index=15;
R(4).index=2;
R(5).index=4;
R(6).index=10;

%plate pairings by region number
if region<=3, 
    P(1)=R(1); 
    P(2)=R(2); 
end;
if region==4, 
    P(1)=R(1); 
    P(2)=R(3); 
end;
if region==5, 
    P(1)=R(3); 
    P(2)=R(4); 
end;
if region==6 | region==7, 
    P(1)=R(1); 
    P(2)=R(4); 
end;
if region>=8 & region<=15, 
    P(1)=R(4); 
    P(2)=R(5);     
end;
if region==16 | region==17,
    P(1)=R(1); 
    P(2)=R(6);     
end;

%Pacific plate is reference
if P(1).index==0; 
    P(1).ws=[0 0 0];
    %P(1).wc=[0 0 0];
else,
    P(1).ws=ws(P(1).index,:); 
    %P(1).wc=wc(P(1).index,:); 
end;
P(2).ws=ws(P(2).index,:); 
%P(2).wc=wc(P(2).index,:); 



%PJH, appears correct 'morvel.m' 
%R.E(2).lon=-112.8; %manual correction to CO lon
